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Abstract. In the present work, we use a photoionisation code to study the H2 emis- 
sion of the Helix nebula (NGC 7293) cometary knots, particularly that produced in the 
interface H + /H° of the knot, where a significant fraction of the H2 1-0 S(l) emission 
seems to be produced. Our results show that the production of molecular hydrogen in 
such region may explain several characteristics of the observed emission, particularly 
the high excitation temperature of the H2 infrared lines. 



1. Introduction 



High-resolution images of the Helix nebula (NGC 7293) have shown that the H2 emis- 
sion arises from its large population of dense globules embedded in the ionised gas (e.g. 
Matsuura et al.|2009 ), the so-called cometary knots (CKs). CKs are structures that re- 
semble comets, particularly in images taken in Ha, [Nil] A65S3, and H2 1-0 S(l) lines. 
The bright cusp points towards the central star and the tail in the opposite direction, 



which can indicate that the excitation is connected with the central star (O'Dell et al. 



2005 2007). 



The H2 emission is intense in a thin layer on the surface of the CKs towards the 



central star. There is no evidence that this emission is produced by shocks (O'Dell 
et al.||2005| |Matsuura et al.]|2007| ). Models of photodissociation regions (PDRs; Cox 
et al.|1998||Q'Dell et al.|200 7 ) are unable to reproduce the high excitation temperature 
of H 2 emission ( -900-1800 K) estim ated by |Cox et al.| ( |1998[ ) and |Matsuura et al. 
(2007). Recently, |Henney et al.| ( [2007] ) showed that advection can cause the ionisation 



and dissociation front to merge, leading to enhanced heating of the molecular gas. 

In this work, we show that the partially ionised region can account for a significant 
part of the observed H 2 emission and naturally explain its high excitation temperature. 
Our models are briefly described in Section 2; a more detailed description will be pub- 
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lished in a forthcoming paper ( Aleman et al.||2010 ). Results are presented in Section 
3. 



2. Models 



We use the photoionisation code Aangaba (Gruenwald & Viegas 1992 1 to simulate 
the ionising spectrum, physical conditions, density of the species, and line emissivities 
around and inside the H + /H° interface of the Helix nebula CKs. The H2 micro-physics 
is included in the code, as described in Alem an & Gruenwald| ( |2004 2010 1. 

We assume that the central star radiates as a blackbody with T+ = 120,000 K and 
L+ = 100 L Q ( |Henry et al.|1999||Q'Dell et al.|2007]> We also assume that the density 



of the diffuse gas is uniform and equal to 50 cm J (Meixner et al. 



abundances for the Helix were determined by Henry et al. ( 1999 1 for He, O, C, N, Ne, 



2005). The elemental 



S, and Ar. For Mg, Si, CI, and Fe, we adopt averages for PNe from Stasiriska & Tylenda 
(|1986|). We use amorphous carbon dust (with 0. l//m radius) for our calculations, but 



as discussed in Aleman & Gruenwald (2004 ) the choice of compound among the ones 
available in the code will not cause significant changes on the results presented below 
nor will affect our conclusions. The distance to the Helix is assumed to be 219 pc 
([Harris et al.|2007). The CKs are simulated as an increase in the density profile of the 



Helix nebula model at a given distance from the central star. The emissivity along the 
radial direction (through the CK symmetry axis) is calculated by the photoionisation 
code. An IDL routine was developed to simulate a three-dimensional CK, allowing the 
calculation of line surface brightness by the integration of the emissivity in the line of 
sight inside a CK, which is assumed to be seen edge on. We construct a grid of CK 
models with different core densities, density profiles, dust-to-gas ratios, and distances 
from the central star. We assume that the density profile has a density increase from the 
diffuse gas to the CK core value within a given distance. We call this region the CK in- 
terface. We study four types of density increase with distance, but here we only present 
the results for the exponential increase (more results will be included in |Aleman et al. 
|2010[ ). Calculation are stopped where the gas temperature, which generally decreases 
with the distance to the central star, reaches 100 K. In each model, the dust-to-gas ratio 
and the chemical composition of the CKs are assumed to be the same as in the diffuse 
gas. 



3. Results 



3.1. Warm H 2 1-0 S(l) emission 

The emissivity of the H2 1-0 S(l) line, as well as for other rovibrational lines, in the 
CKs is important in a warm region, where temperatures are between 300 and 7000 K. 
In the region considered in this work, the peak in the 1-0 S(l) emissivity occurs where 
the density is around 40% of the core density. The contribution of colder regions should 
be more important for pure rotational lines of the v = level. This component of the 



H2 emission may explain the excitation temperatures around 900-1800 K found by [Cox 
etaLl ( fl99"8] ) and |Matsuura et al.| ( [2007l ). 

The left panel in Fig. T shows the H2 excitation diagram. Observational values 



were obtained by Matsuura et al. (2007 1. Values calculated with an appropriate model 



are also included. The agreement between the excitation temperatures of the model and 
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Figure 1. Left Panel: H2 excitation diagram. The effective column density was 
calculated where the 1-0 S(l) surface brightness is maximum. Open symbols repre- 
sents observations and filled symbols models. Lines are Boltzmann distributions for 
the temperature indicated. Right panel: H2 1-0 S(l) surface brightness of a cometary 
knot as a function of its distance to the central star. Sets of solid, dashed, and dot- 
dashed curves represent models with interface thickness of 0.5, 0.2, and 0.01", re- 
spectively. Different curves within each set represent CK radius of 0.5, 1 .0, and 2.0", 
with the surface brightness increasing for larger CK radius. Dots represent measured 
values. The error in the distance is estimated assuming that the Helix symmetry axis 
is inclined 37 with respect to the line of sight 



the observations is evident. Lines represent Boltzmann distributions for three different 
temperatures as indicated within the plot. The column densities obtained from the lines 
of the bands 1-0 and 2-1 are well represented by a excitation temperature of approxi- 



mately 2000 K. A similar value was obtained by Mat suura et al.| ( |2007| ). The column 
densities obtained from lines 0-0 S(2) to S(7) are well represented by Boltzmann distri- 
bution at a temperature of 1000K, which is close to the excitation temperature of 900K 



obtained by |Cox et al. (1998 ) from ISO observations of the Helix. 



3.2. H 2 1-0 S(l) surface brightness 

The right panel of Fig. [T] shows the H2 1-0 S(l) line surface brightness as a function 
of the CK distance to the central star. Dots represent measurements for some repre- 
sentative CKs of the Helix nebula. We identified 10 isolated CKs commonly detected 
in Ha and H2 images. We measured 2.12 /mi H2 intensities from the images obtained 



by |Matsuura et al. (2009). To calibrate the data for point sources on the local scale, 
we use five stars within the observed field to measure the zero-pint. We assume that 
the 2MASS ^'-magnitude of these stars are the same as the magnitudes in H2 -filter. 
We apply the 25-pixel radius of aperture photometry and take the 35-50 pixel ring as 
background measurements. The pixel scale is 0.117". 

Curves in the right panel of Fig. [T]show the H2 1-0 S(l) surface brightness as a 
function of the distance from the central star for CKs with different interface thickness 
and radius. The surface brightness was calculated was averaged over the same aperture 
as the measurements to allow direct comparison. The surface brightness decreases with 
the decrease of interface thickness, the increase of the CK radius, and the increase of the 
distance from the central star. The observed surface brightness also presents a decrease 
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with distance to the central star trend. The interface can account for the whole or a 
significant part of the observed surface brightness. 

An important parameter to the ionization structure of the CKs is the distance from 
the central star, since the ionising spectrum may change significantly with the optical 
depth. CKs farther from the central star have smaller ionised zones. If the CK is beyond 
the Helix ionisation front, there is practically no ionised region and the intensity of 1- 
S(l) drops dramatically, since there is not enough radiation or temperature to excite 
significantly the upper vibrational levels of the molecule. Our results support that the 
central star's radiation field plays a major role in the CKs H2 emission. 

We also study the effect of the dust-to-gas ratio and hk . The 1-0 S(l) peak bright- 
ness is slightly higher for models with higher dust-to-gas ratio (the increase caused 
by changing the dust-to-gas ratio from 10 -3 to 10~ 2 is about 20%) and higher hk (the 
difference between models with hk = 10 5 and 10 6 cm" 3 is up to 40%). 



As pointed out by Burkert & O'Dell ( 1998 1, the interface between the diffuse gas 



and the CK core may provide important clues about the mechanisms that shapes and 
sustain the CKs. Our models show that there are significant differences in the results 
depending on the assumed density profiles of this region. Images that could resolve the 
interface of the CK in great detail are then essential to improve or knowledge about the 
CKs and PNe. 
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